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The efficiency of frequency up and down conversion of light in an atomic ensemble, with a diamond 
level configuration, is analyzed theoretically. The conditions of pump field intensities and detunings 
required to maximize the conversion as a function of optical thickness of the ensemble are determined. 
The influence of the probe pulse duration on the conversion efficiency is investigated by numeric 
solution of the Maxwell-Bloch equations. 
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I. INTRODUCTION 

The frequency conversion of light fields has been an 
important theme in optical physics for around half a cen- 
tury. In quantum information physics the conversion of 
single photons to and from the telecom wavelength band 
is a topic of more recent vintage, and is motivated by 
the desire to minimize optical fiber transmission losses 
when distributing entangled states over distant quantum 
memory elements in a quantum repeater 

An associated technical problem is that telecom light 
is not readily stored in ground level atomic memory co- 
herences. Retrieval processes in atomic ensembles , for 
example using electromagnetically induced transparency 
or more specifically the dark-polariton mechanism 
generate shorter wavelength radiation correlated 
to the stored atomic excitation by Raman scattering. 
Such radiation, optically resonant to the ground level 
of typical atoms and ions, has been retrieved in numer- 
ous experiments [Bl-tlil] . An important advance would in- 
volve generation of atomic memory coherences quantum- 
correlated with telecom wavelength radiation, thereby 
minimizing transmission losses over long distances. Re- 
cently there has been a breakthrough in this direction 
using a pair of cold, nondegenerate rubidium gas sam- 
ples [15|. The stored excitation is correlated with an 
infrared field (idler) in one gas sample, and the idler is 
then frequency converted to a telecom wavelength signal 
field in the other ensemble. The frequency conversion 
mechanism involves the diamond configuration of atomic 
levels shown in Fig. 1. 

In a probabilistic protocol it is important to maximize 
all efficiencies, e.g., fiber transmission single-photon de- 
tection and quantum memory lifetime [16]. In the present 
work we investigate the efficiency of frequency up- and 
down- conversion in the diamond atomic configuration 
[l7l [l9j . as a function of the ensemble's optical thick- 
ness, and the intensity and detuning of the pump fields 
involved in the near-resonant, four-wave mixing process. 

The remainder of this paper is organized as follows. 
In Sec. II we discuss the four wave mixing process and 
present solutions for the up- and down- converted fields; 
the dressed state picture is used as a guide to under- 
stand the characteristic features of the absorption and 



signal-idler field coupling, as a function of the micro- 
scopic interaction parameters. In Sec. Ill we present 
results of an optimization of the conversion efficiency as 
a function of the optical depth of the atomic ensemble. 
In Sec. IV we investigate the effects of finite pulse dura- 
tion by integrating the Maxwell-Bloch equations for the 
system. Section V presents our conclusions. The deriva- 
tions of the Maxwell-Bloch and parametric equations are 
relegated to the Appendix. 



II. UP AND DOWN CONVERSION 
EFFICIENCY 

In this paper, we consider a cold and cigar-shaped 87 Rb 
atomic ensemble with co-propagating light fields similar 
to the experimental setup in Ref. [15| . 

The conversion scheme shown in Fig. Q] involves two 
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FIG. 1. The diamond configuration of atomic system for con- 
version scheme. Two pump lasers (double line) with Rabi 
frequencies fi a ,fij, and propagated probe fields (single line) 
Ef,Ef interact with the atomic medium. Various detun- 
ings are defined in the Appendix and the atomic levels used 
in the experiment QJ] are (|0), |1), |2), |3» = (|5S 1/2 ,F = 
1), |5P 3/2 , F = 2), |6Si/ 2 , F = 1), |5P 1/2 , F = 2)). 
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pump lasers with frequencies uj a and uib, respectively; 
their Rabi frequencies are given by VL a and f^. Two 
weak probe fields, signal and idler, with frequency uj s and 
uji, respectively, propagate through the optically thick 
atomic medium. Unlike the cascade driving scheme, 
where two-photon excitation generates a photon pair 
spontaneously 17], pump laser b experiences a trans- 
parent medium if both the signal and idler fields are in 
the vacuum state. With an incident signal held, four 
wave mixing with the pumps generates an up-converted 
idler field, while an incident idler field generates a down- 
converted signal. 

The Maxwell-Bloch equations for the interacting sys- 
tem of light and four light fields is derived in the Ap- 
pendix. By linearizing the equations with respect to 
the signal and idler field amplitudes, and adiabatically 
eliminating the atoms, one arrives at coupled paramet- 
ric equations for the signal and idler fields. We discuss 
their solution in this section, and leave numerical solu- 
tions of the Maxwell Bloch equations to Sec. IV. The cal- 
culation of conversion efficiencies can also be carried out 
with the quantized Heisenberg-Langevin version of the 
coupled parametric equations. The resulting conversion 
efficiencies are identical to the semi-classical treatment; 
the additional quantum noise contributions vanish as the 
1 2) 1 3) transition driven by pump laser b has vanishing 
populations and atomic coherence. A similar simplifica- 
tion occurs in the calculation of the storage efficiency of 
spin waves in a system of atoms in the A configuration 



The co-moving propagation equation for c-number sig- 
nal and idler fields (respectively, Ef and Ef) under en- 
ergy conservation (Alo = u a + lu s — uib — uji = 0) and 
four- wave mixing conditions (Afc = k a — k s + kb — ki = 0) 
are 
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For up-conversion, (Ef (0) = 0) symmetry gives the cor- 
responding coefficients 
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The up and down conversion efficiencies differ only 
by the interchange ft, ■<-> k s . In the strong parametric 
coupling regime where |«i|,|At s | S> |ci!.j|, |/3 S |, the coeffi- 
cients can be simplified to ?7 U — smh(^/ K s KiL) , r/d ~ 

smh(y/K s KiL) and T u = Td ~ cosh(^/ n s KiL). 

Under the further assumptions a, = /3 S — and Ki , k s 
are pure imaginary, we find ?7 U = ?/d = sin [Im(«; s .L)] 
and T u = — cos 2 [Im(K s L)]. This result was re- 
cently derived by Gogyan using a dressed-state approach 
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The coupled equations are similar to those found for the 
double A system [2(3, [13 ■ The self-coupling coefficients 
/3 S , ai and parametric coefficients defined in 

the Appendix. The solution, under conditions of down 
conversion with boundary condition Ef(0) = 0, is 
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where w = \/ q 2 + and q = (— a,; + /3 s )/2. We 

define the down conversion efficiency ?yd and transmission 
of input idler field Td by 




FIG. 2. Dressed-state picture from the perspective of the 
probe idler transition between atomic levels |0) and |3). Two 
strong fields Q a ,flh shift the levels with energy AE a ,t and 
wavy lines represent the idler field resonances. 
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J, in the case of resonant pump fields Ai — At, = 
!] . In this ideal limit there is a conservation condition 
Vu + T n = r)d + Td = 1- The parametric coupling coeffi- 
cients are not identical, but in the regime of strong cou- 
pling they approach each other. As noted by Gogyan, 
when the pump-a intensity is large (fl a >> | Ai 1,703) 
the |0) f-> 1 1) is saturated, the atomic coherence is neg- 
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tively, in the limit fib >> fla, ^32 and |Ai| >> 701 the 
atomic coherence of |0) O |1) dominates and once again 
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The ac-Stark splitting induced by the pump lasers 
shifts the resonant absorption condition for the idler 
and signal fields. The idler and signal experience 
resonant absorption at the transition frequency of 
the dressed atom. The corresponding transitions for 
the idler are shown in FigJ5] The bare states are 

Ariv/Af 



shifted by AE a 



/2 and AE b 



A b ± VA 



ml 



/2, respectively. Note that our Rabi 

frequencies are smaller by a factor 2 than the standard 
definitions to avoid a plethora of prefactors in the equa- 
tions of the Appendix. 



For resonant pump fields, AE a ^ — ±fl a ,b . The 
idler transition resonances are at Auii — — (fl a + fib), 
— \fl a — fl b \ ! |n a — , (fl a + fib) and these delineate 
three windows separated by these four absorption peaks. 
For fl a > fif, the centers of these windows are at — fi , 
0, and fl a , respectively. Choosing the idler detuning 
Aoji = ±fl a as in [23j . the idler interacts with the atomic 
medium at the center of the left or right window. 



As an example of the strong coupling windows cre- 
ated by intense pump lasers, we show in Fig. [3] the self- 
and cross-coupling coefficients for the signal and idler 
fields as a function of the idler frequency. Note that the 
corresponding frequency of signal field is determined by 
Auj s = Aui— Ai+At,. The dimensionless quantities onL, 
(3 S L, and n s L are shown under the conditions of maxi- 
mum conversion efficiency to be discussed in the next 
section. We choose the optical depth (opd) paL =150 
where p is the number density, a = 3A 2 /(47r) the reso- 
nant absorption cross section, and L the atomic ensemble 
length in the propagation direction. Three parametric 
coupling windows are separated by two strong absorp- 
tion peaks on the left and two relatively weak ones on the 
right. The imaginary part of the self-coupling coefficients 
are seen to vanish in each window at a certain point, 
while the real parts are small away from resonances. At 
the same time the cross-coupling coefficients have a large 
imaginary part. The positive gradient of Im(/3 S L) and 
Im(aiL), inside the windows is indicative of normal dis- 
persion. 



III. OPTIMAL CONVERSION EFFICIENCY 

It is important to ascertain the parameters which al- 
low maximum efficiency of conversion due its potential in 
practical quantum information processing. In principle 
we need to search the three parametric coupling windows 
to find the optimum conditions for an atomic ensemble 
of a given optical thickness. 

In the previous section we have discussed how three 
parametric coupling windows appear for some particular 
values of pump laser parameters. In the search for the 
maximal conversion efficiency, five parameters fl a , fib, 
Ai, Ab, and Au>i are varied to maximize the conversion 
efficiency for a fixed optical depth of atomic ensemble, 
using functional optimization. 

The optical depth poL appears through the depen- 
dence on atomic number N in the Arecchi-Courtens co- 
operation time T c [25j 
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In Fig. 21 we show the maximum of down conversion effi- 
ciency using Eq. Q for different optical depths from I to 
300. The maximum is found by varying five parameters 
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FIG. 3. (Color online) Self-coupling coefficients /3 s ,ai and 
cross-coupling coefficient k 3 . Dimensionless quantities (a) 
/3 S L, (b) ctiL and (c) k s L with real (solid blue) and imagi- 
nary (dashed red) parts are plotted as a dependence of idler 
detuning Auii [same label in (b)] showing four absorption 
peaks to construct three parametric coupling windows. A 
black dashed-dot line of the constant 7r/2 is added in (c) to 
demonstrate the crossover with Im(K s L) indicating the ideal 
conversion efficiency condition in the left window. The pa- 
rameters we use are ($l a , fib, Ai, A5, Auii) = (33, 20, 39, 2, 
—21)703 for optical depth paL = 150 with L — 6mm. Vari- 
ous natural decay rates are 703 = 1/27. 7ns, 701 = 1/26. 24ns, 
712 = 703 /2.76 and 732 = 703/5.38 
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mentioned above and the conversion efficiency reaches 
100% asymptotically when the optical depth becomes 
larger. In the strong parametric coupling regime as 
we discussed in the previous section, rjd ~ sin 2 [Im(K s L)] 
and it has a maximum when Im(«; s L) = -|, see Fig|3] 
Since Im(«; s L) is proportional to optical depth and in- 
versely proportional to the Rabi frequencies of the driv- 
ing lasers, an order of magnitude estimate of the opti- 
cal depth necessary for near unit conversion efficiency is 
opd- § fVfc/703 » L 

The behavior of the cross-coupling coefficient 3(k s L) 
as a function of idler detuning indicates where large con- 
version is to be found, as a comparison with Fig(5]shows. 
The maximum efficiency of about 0.92 is located in the 
left parametric coupling window at the intersection of 
3(k s L) and J. Inside the windows the trade-off be- 
tween conversion and transmission is clear. In the re- 
gion where absorption is large, on the sides of the win- 
dow (especially for the left window), the efficiency and 
the transmission are both low although the valley in con- 
version efficiency corresponds to a peak in transmission 
as expected in parametric coupling. The transmission 
approaches unity when the incident idler field is far off- 
resonance. 

We note that the symmetry (Ai, At, Aw,) — > 
— (Ai, Aft, Aw,) gives degenerate optimal conversion con- 
ditions. 



IV. PULSE CONVERSION: SOLUTION OF THE 
MAXWELL-BLOCH EQUATIONS 

The effect of finite-duration input probe pulses, which 
are often employed in practice, can be assessed by nu- 
merically solving the Maxwell-Bloch equations for the 
coupled atoms-fields system. The characteristic scales of 
time and length are given by the Arecchi-Courtens time 
T c and L c — cT c , respectively, which are inversely pro- 




100 200 300 400 500 600 

optical depth 

FIG. 4. (Color online) Down conversion efficiency rjd vs op- 
tical depth (opd) from 1 to 300. Each dotted point is the 
maximum for five variational parameters Q a ,Qb, Ai, At, and 




-80 -60 -40 -20 20 40 60 80 



AB Y%3 



FIG. 5. (Color online) Conversion efficiency r)d, f] u and trans- 
mission Td vs Au)i for opd=150. 7/d and rj u are indistinguish- 
able and shown in solid red line and Td is in dashed blue 
line. High transmission efficiency corresponds to low conver- 
sion efficiency indicating the approximate conservation condi- 
tion within each parametric coupling window. The maximum 
conversion efficiency is found in the left window at around 
Acji = — 2O703 and other relevant parameters are the same as 
in FigO 



portional to the square root of the opd. The cooperative 
electric field is the product of the atomic number and the 
idler electric field per photon, i.e., E c = ^ pHuii / (2eo) ■ 

Scaling the space, time, and electric field amplitude ac- 
cordingly, indicated by tildes, the light propagation equa- 
tion becomes 
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A consistent scaling can also be applied to the atomic 
dynamical equations presented in the Appendix. 

The Maxwell-Bloch equations were integrated with 
a semi- implicit finite difference method 26]. The mid- 
point integration method is stable and has high accuracy 
without sacrificing memory for finer grids 27] . The algo- 
rithm has been tested by comparing with the parametric 
equations' solutions in appropriate limits, and these solu- 
tions are recovered when fine enough grids are employed. 

To illustrate the influence of finite pump pulse dura- 
tion, we compute the down conversion efficiency 



f\E+(z = L,t)\*dt 
f\E+(z = 0,t)\*df 



(8) 



In Fig. |51 we show the computed values of r]d for 
two different input idler pulse durations. We fix the 
opd= 150 and use the near optimum parameters (fi a , fib, 
Ai, A 6 , Aw,) = (33 , 20, 39 , 2, -21)7 03 determined from 
the coupled parametric equations. The temporal shape of 
the pump laser intensities are also shown. Pump laser b 
is taken to be continuous wave, while pump a is a square 
pulse with duration large enough to completely overlap 
the input idler pulse. To compare with the steady state 
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FIG. 6. (Color online) Time-varying pump fields of 
Rabi frequencies Q a ,b(t) and down-converted signal intensity 
[\E+(t,z = L)\ 2 ] from an input idler pulse [\E+(t,z = 0)| 2 ]. 
Pump-b (dotted green) is a continuous wave and pump-a 
(dashed black) is a square pulse long enough to enclose input 
idler with pulse duration (a) 100 ns and (b) 15 ns (dashed- 
dot blue). Output signal intensity (solid red) at the end of 
atomic ensemble z = L is oscillatory due to the pump fields. 
The square pulse in rising region (t r — \ < £ < ir + 4f) has 
the form of |[1 + sin( 7r ft~* r ? )] that in (a) (U,t s ) = (10, 10)ns 
for pump-a and (t r ,t s ) = (20, 20)ns for input idler; (b) 
(t r ,t a ) — (10, 5)ns for pump-a and (t r ,t s ) = (15, 10)ns for 
input idler where t r is the rising time indicating the center of 
rising period t s . Note that the falling region of square pulse 
is symmetric to the rising one. 



V. CONCLUSION AND DISCUSSION 

We have studied light frequency conversion in an 
atomic ensemble with a diamond configuration of atomic 
levels such as 87 Rb. The motivation stems from the need 
to efficiently convert light resonant with ground state 
transitions (storable in the sense of quantum memories) 
to and from the telecom wavelength band for low-loss 
quantum network communication. The optically thick 
atomic sample is driven by two strong co-propagating 
pump fields, and a probe idler or signal field depending on 
whether we consider down- or up-conversion. Paramet- 
ric equations for the probe fields are derived and used to 
compute conversion efficiencies. By performing a global 
search we find conditions of pump Rabi frequencies, de- 
tunings and signal / idler input frequency to maximize the 
conversion efficiency as a function of optical depth of 
the ensemble. Only in the limit of very large optical 
depth does the maximum efficiency approach the ideal 
strong coupling result (23[. Under conditions routinely 
obtained in cold, non-degenerate rubidium gas, with opd 
~ 100—200, optimal conversion efficiencies of the order 80 
to 90% are predicted. Numerical solution of the Maxwell- 
Bloch equations confirms the solution of the parametric 
equations in the limit of long pulse duration, and indi- 
cates that for shorter pulses, pump pulse induced modu- 
lation may reduce the conversion efficiency. 
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solutions, we choose the Rabi frequency of idler as O.I703 
which is small compare to those of the pumps. We find 
that the conversion efficiency is reduced for shorter idler 
pulse inputs. A 100 ns idler pulse is long enough that it 
has a almost the same maximum conversion efficiency of 
0.92 as in Fig|4]for opd=150. While for the shorter idler 
pulse of 15 ns, the signal develops significant temporal 
modulation and this reduces the conversion efficiency, al- 
though it is still quite appreciable. The modulation fre- 
q uency is at the generalized Rabi frequency of pump-a 
\J A\ + 40^. We note the characteristic time and space 
scales of the calculations are T c = 0.086 ns and L c = 26 
mm for a moderate atomic density p = 1.7 x 10 11 cm -3 
and L = 6 mm. 

The grid size for dimensionless time At = 0.5 and space 
Az = 0.001 were chosen for both 100 and 15 ns idler 
pulse durations and the convergence is reached with an 
estimated relative error less than 1%. 



Appendix A: Maxwell-Bloch and parametric 
coupling equations 

To derive the coupled Maxwell-Bloch equations it is 
convenient to employ a quantized description of the elec- 
tromagnetic field [28j and use Heisenberg-Langevin equa- 
tion methods, and then invoke a standard semiclassi- 
cal factorization assumption. The propagation length 
L is discretized into 2M + 1 elements. The posi- 
tive frequency component of the electric field operator 

is given by E+(z) = J2n=-M \f^^e^ ks+k ^ z c n where 

[cVj, C n f] — Snn f , k n — ^ 1 ^s.n — ^ s ~\~ k n C , 71 — 

—M, M and uj s = k s c is the central frequency. Define 
the local boson operators a; = ^ 2 l I+1 Yln=-M c n e lknZl 

where [a;,a|,] = 6u> . Similar definitions hold for the sig- 
nal, s, and idler field, i, which carry an additional index 
in the following. 

The Hamiltonian for the interacting system depicted in 
Fig Q] is given by, (we ignore the interactions responsible 
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for atomic spontaneous emission for the moment) 

H = H + Hi, (Al) 

where 
H 
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where at, 



the Rabi frequency 



(&)(*) = fa,(b){t)d w ,(23)£(k a ^ b) )/(2h) is half the stan- 
dard definition and f a ,(b) is a slowly varying tempo- 
ral profile without spatial dependence (ensemble size 
much less than pulse length). The dipole matrix 
element d m n = (m\d\ri), coupling strength g s ,(i) = 

d2i,(30)£( k s.(i))/h, £(k) = ^/hw/2e V and z p = 
pL/(2M + 1), p = -M, M. The matrix uj w = 
E n =-M« k n e lk ^ z '- z ''^/(2M+l) accounts for field prop- 
agation by coupling the local mode operators. 



The dynamical equations including dissipation due 
to spontaneous emission may be treated by standard 
Langevin-Heisenberg equation methods [29j], and we de- 
fine 7y as the natural transition rate from \j) — > Since 
we are interested in a semiclassical description, we replace 
the field operators by c-numbers in the Langevin equa- 
tions, and drop the zero-mean Langevin noise sources. 
All atomic spin operators are also replaced by their ex- 
pectation values. Finally, in the co-moving frame coor- 
dinates z and t — t — z/c the atomic equations are 
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where 72 = 712 + 732, P = e lAkz lAwt ) the four- wave varying atomic operators are defined 
mixing mismatch wavevector Afc = k a — k s + k b — fcj, the 

frequency mismatch Aw = w a + w s - uj b - uji = Ai - _ J_ , e _ lfeaZi+JWat _ J_ ijfe.^+w.t 
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where N Z (2M + 1) = TV . 
The field equations are 
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where the field operators are defined as 
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For energy and momentum conservation (P = 1), and 
in the weak field limit, we solve the atomic equations in 
steady state after linearizing with respect to the probe 
fields 
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expressed in first order of fields and &oi, O32 m zeroth 
order. For population operators, we solve them in the 
zeroth order of fields and the nonzero steady states of 
population and coherence operator are (s denotes steady 
state solution) 
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Substitute the latter results into Eq. (|A8[) . and solve for 
CT12 and <7o3- The parametric coupling equations for the 
signal and idler fields become 



where 
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^-Et = K i E+ + a i Ef 
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